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Abstract 

In this paper we carry out a stability analysis of the Bloch states of a two-component Bose- 
Einstein condensate confined to a Id optical lattice. We consider two concrete systems: a mixture 
of two hyperfine states of Rubidium-87 and a mixture of Sodium-23 and Rubidium-87. The former 
is seen to exhibit similar phenomena to a single component condensate while the latter also suffers 
an instability to phase separation at small Bloch wave vectors. It is shown that sufficiently deep 
optical lattices can remove this latter instability potentially allowing imiscible cold atoms species 
to be held in intimate contact and transported within an experimental system. 
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I. INTRODUCTION 



Bloch states of weakly interacting Bose condensates have been the subject of recent 
experimental and theoretical study (see for a review). While these states are, in many 
respects, similar to those of electrons in solids, they differ in one crucial respect - they are not 
necessarily stable . Isolated, weakly interacting atomic condensates are well described 
by mean field theory in the form of the time-dependent Gross-Pitaevskii (GP) equation. In 
the presence of an optical lattice this has stationary Bloch solutions of the form 

tf(x,t) = e-^ t/h e ikx f{x) (1) 

where / (x) has the periodicity of the lattice and \x is the chemical potential for the atoms. 
Many of the phenomena associated with these Bloch states are familiar from the theory of 
electrons in solids (band structure etc.). Indeed, cold atoms systems readily exhibit phe- 
nomena which are hard to observe in electronic systems such as Bloch oscillations. However, 
unlike electronic states in a crystalline lattice, the Bloch solutions of the GP equation are not 
necessarily stable. As shown by Wu and Niu and Machholm, Pethick and SmithQ,!^, there 
are two types of instability. Firstly there is an energetic instability associated with Bloch 
states which are not local minima of the mean field energy, secondly there is a dynamical 
instability associated with the exponential growth of small perturbations around the Bloch 
state. The former instability should not be visible in systems which are perfectly described 
by the GP equation, which conserves the mean field energy, but will have an effect once 
dissipative processes become active. The dynamical instability is always significant unless 
the time for which the system is described by the unstable Bloch state is much shorter than 
the growth time for the most unstable mode. 

In this paper we consider the dynamical stability of a system consisting of two distinct 
atomic species. The laser generating the standing wave optical lattice is supposed to have 
been chosen so that it is blue detuned from the nearest resonance of one species but red- 
detuned from the nearest resonance of the other. This means that the two species will see 
opposite potentials as one is attracted to the nodes of the standing wave while the other 
is attracted to the antinodes. We also suppose that the system is strongly confined by a 
cylindrically symmetric magnetic trap to ensure one dimensional behaviour. 

In section 2 we review the mean field theory for such a system and the 2-component 
GP equation. In section 3 we consider the Bloch states of the 2-component condensate. In 



section 4 we review the linear stability analysis of the mean field theory. In section 5 we 
consider two specific systems: firstly we consider the two species to be 2 hyperfine states of 
Rubidium; next we consider a mixture of Rubidium and Sodium atoms. We will see that 
these two examples exhibit rather different behaviour, the former is qualitatively similar to 
a single component system while the later exhibits quite different behaviour. In section 6 
we will discuss the results obtained. 



II. MEAN FIELD THEORY FOR 2-COMPONENT SYSTEM 



The mean field theory for a two component system is well described in the book by 



Pitaevskii and Stringari 



The main result is that the two component condensate wave 



function, tyj (x, t) satisfies the coupled GP equations 
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where rrij is the mass of an atom of species j and Vj is the optical lattice potential seen 
by species j. The nonlinear terms arise from the treatment of the atomic collisions within 
mean field theory. For the low temperatures relevant to ultra-cold atom experiments it is 
sufficient to treat the atoms as point scatterers with strengths given in terms of the relevant 
s-wave scattering length within the lowest Born approximation. This gives 
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for 2-body processes in which like atoms scatter and 
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for 2-body processes in which unlike atoms scatter. The t^'s are form-factors accounting 
for the finite extent of the wavefunctions transverse to the optical lattice. For a harmonic 
radial trap of the form U (r) = C 2 r 2 /2 we have 

J " = WW) (6) 

where the radial oscillator lengths for the two species are I 2 = h/Cyfrnl. The expression 
for g>i2 involves the reduced mass m\i = (rrii 1 + m^j relevant for the "unlike" 2-body 
problem. 



Stationary solutions of these equations of the form tyj (x, t) = e w*/ n ^- (#) satisfy the 
time-independent GP equations 

h 2 d 2 ^ 
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in which the chemical potentials of the components appear as the solution to a non-linear 
eigenvalue problem. 

The time-independent GP equation can be obtained variationally from the mean field 
energy functional 
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where, as usual, the chemical potentials arise as a Lagrange multipliers enforcing the nor- 
malization conditions 

Jdx\^\ 2 = N 3 (11) 

for a system with Nj atoms of species j. 

The chemical potential of a stationary state is related to the mean field energy via the 
relation 

£ = (fuN, + fi 2 N 2 ) - \gij J dx \A\ 2 l^-l 2 • (12) 

ij 

We consider a system in an infinite optical lattice with 

Vj (x) = Wj cos (kx) (13) 

where k = 2n/d where d — A/2 is the period of the lattice (half of the wavelength of the 
laser generating the standing wave). 

III. BLOCH STATES 

Next we seek Bloch states of the form 

0f (x) = e ik *f 3 (x) (14) 



4 



where 

f j (x + d) = f j (x) . (15) 

We employ the same basic method as Machholm and Smith [3] , adapted to the case of two- 
components. If rij is the number of atoms of species j per unit length then the normalization 
condition on fj(x) becomes 

1 rd/2 

dx \fj(x)\ = rij . (16) 



d i-d/2 

We also define the energy per unit length 



E 
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The periodicity of fj (x) allows us to write 

./}(•'') (19) 



so that 



and 
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As found by Machholm and Smith for the single component case, we can restrict our atten- 
tion to real values of the / parameters and truncate summations over the Fourier label to 
values less than a cut-off S. We then minimize £ as a function of the 2 x (2S + 1) / pa- 
rameters subject to the constraint (}20|) to find the Bloch wave functions for each component 
for a grid of points in the range — k/2 < k < k/2 (i.e. within the first Brillouin zone of the 
optical lattice). 



IV. STABILITY ANALYSIS 



Given the Bloch state for a particular value of k, we can ask whether it is stable with 
respect to small fluctuations. We assume an initially small, generic fluctuation, so that the 
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condensate wave function has the form 



*j (x, t) = e-^ t/n (fa (x) + (x, t)) 



(23) 



which we substitute into the time-dependent GP equation. Dropping terms non-linear in 
the fluctuations then gives 
88%! 
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As expected, the <5\l/'s are coupled to their complex conjugates and we must decouple them 
using a classical version of the Bogoliubov transformation, as developed by Pitaevskiij^]. 
We set 

6*j (x, t) = e ikx (e^-^Uj (x) + e-^-^v* (3)) (28) 

where Uj (x) and Vj (x) all have the periodicity of the lattice. Substituting into the time 
dependent equation gives, after some manipulation, an eigenvalue problem for the u(x) and 
v(x) functions of the form 
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Substitution of the truncated Fourier expansions 
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allows this to be turned into an 4 (2v + 1) x 4 (2v + 1) generalized matrix eigenvalue problem. 
The resultant eigenvalue problem is non-hermitian and hence need not, in general, have real 
eigenvalues. As shown in fj the eigenvalues either come in real pairs ±u or in sets of 
four complex eigenvalues ±u/ ± iu" . As usual in the Bogoliubov method these are not all 
independent. For real eigenvalues only modes with positive frequency need be considered. 
For the case of complex eigenvalues we may discard those with negative real parts. The 
presence of an imaginary part to the eigenvalue indicates an instability of the underlying 
Bloch state, since there will be a mode which, at least initially, grows exponentially in time 
at the rate to" . 

The eigenvalues for each q value form a set corresponding to the Bloch bands of the 
linearized fluctuation problem. Only the lowest two fluctuation bands exhibit non-zero 
imaginary parts to the eigenvalues, so we focus attention on these. We define the instability 
of the system with respect to modes with wave-vector q as 

9(q)= sup {u'j(q)\ (35) 

V. RESULTS OF STABILITY ANALYSIS 

Here we will present results for the stability of the Bloch states of 2 systems. We will 
compare and contrast the two systems in the next section. 

Firstly, we consider a system in which the two hyperfine states \F = l,rriF = —1) and 
\F — 2, mp = 2) of Rubidium-87 are cooled in a magnetic trap to form a condensate. Such a 

n 

two-component condensate was first prepared experimentally by Myatt et al[7J. We suppose 
that the atoms are confined by a strong axial magnetic trap and a longitudinal optical lattice. 
We further suppose that the laser is tuned between resonances of the two species such that 
Vq = Wi = —w 2 and we choose a typical value for the depth of the lattice V = 0.765-Er 
where 

E "=2^ < 36 > 
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FIG. 1: The Bloch wavefunctions fj(x) at k = (dotted) and k = k/2 (full) plotted against kx 

is the recoil energy associated with a Rb atom absorbing a photon from the laser generating 
the standing wave. We also assume equal densities of the two species, n\ = n<i = n and 



so that, apart from the opposite signs of the v/s, the two species are virtually identical. In 
figure show the form of the Bloch functions fx and fi for both k = and k = k/2 (zone 
boundary). As can be seen, the k = Bloch states are lightly modulated by the optical 
lattice with one species attracted to the nodes of the optical standing wave and the other 
attracted to the antinodes. The k = k/2 Bloch states show much stronger modulation with 
the wavefunctions of the two species vanishing at the nodes and antinodes respectively. 

Figures 121 and El show colour maps of the imaginary part of the fluctuation frequencies, 
uj'{(q) and oo'^q) respectively, for the lowest two fluctuation bands as functions of the Bloch 
wave-vector k and the fluctuation wavevector q. As can be seen, both of the lowest bands 
exhibit behaviour similar to that of a single component condensate. The Bloch states are 
stable for k < k/4. States with higher Bloch wave-vectors are unstable with respect to 
fluctuation modes with wave- vectors around q = k/2: effectively a period doubling modu- 
lational instability. As k is increased, the unstable modes move to longer wavelengths with 
a greater range of g's being unstable. All of this behaviour is qualitatively similar to the 
behaviour of a single component condensate. 



ngi2 = Q.099E R 



ng n = ng 22 = 0-lE R 



(37) 
(38) 
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FIG. 2: The imaginary part of the fluctuation mode frequency, $suj\ (q) for < k < k/2 and 
< q < k/2 (colour online). 




FIG. 3: The imaginary part of the fluctuation mode frequency, Qu! 2 (<?) f° r < A; < k/2 and 
< q < k/2 (colour online). 

Now we consider a different system: a mixture of Sodium and Rubidium atoms, similarly 
confined to a Id optical lattice tuned between resonances of the two species. We take 
typical values Vo = w\ = —w 2 = 0.603i?R and again assume equal densities n\ = n<i = n 
such that ng n = 0.047E R , ng 2 2 = 0.019£ , r and ngi 2 = 0.048-Ej?. Figure |U shows the Bloch 
wavefunctions for this system at k = and k = k/2. This system is much less symmetric 
and the wave functions for species 1 are more spread out due to the larger intra-species 
interaction. Once again the Bloch state at the zone boundary, k = k/2, has nodes for both 
species. FigureElis a colour map showing, for each Bloch wavevector, k, and each fluctuation 
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FIG. 4: The Bloch wavefunctions fi(x) and ^(s) for a Rb-Na mixture for k = (dotted) and 
A; = k/2 (full). 
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FIG. 5: Instability 6{q) of a Rb-Na condensate as a function of k and g for v = 0.6085Er (Colour 
online) . 

wavevector, q, the instability 0(q). As can be seen, this map is qualitatively different to that 
of a single component system. The map has two regions, small k and larger k. The large 
k behaviour is similar to the one component case in that instability sets in at k w k/4 for 
modes with q = k/2. At larger k the dominant (i.e. most rapidly growing) mode moves 
to longer wavelengths. Unlike the single component case, the short wavelength modes do 
not become stable as k increases and at k — k/2 all g's are unstable. The small k regime 
is quite unlike the single component systems, there are long wavelength unstable modes 
even at k — which persist up to a maximum k, with the dominant fluctuations moving 
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towards shorter wavelengths. This leaves a narrow window of Bloch wave vectors for which 
the system is stable. 

The origin of the instability at small k lies in the fact that this system is unstable even in 
the absence of an optical lattice exhibiting a strong tendency to phase separate into single 
component domains. As shown by 0, 0] the condition for the mode with wavevector q to 
be stable in the absence of a lattice is 



so that if #11(722 < 9i2 only modes with sufficiently large q are dynamically stable. It is clear 
that the parameters for the Na-Rb mixture do not satisfy this stability criterion. 

VI. USING AN OPTICAL LATTICE TO STABILIZE TWO-COMPONENT CON- 
DENSATES 

The depth of the lattice potential used above was chosen rather arbitrarily. We expect 
that increasing the depth of the optical lattice should enhance the stability of the NaRb 
system at low k because it confines the two components in different places. In particular, a 
very large value of Vo should lead to an array of pure phase domains with the same periodicity 
as the optical lattice. In order to see how deep an optical lattice is required to stabilize the 
system we have carried out the stability analysis for the k = Bloch states for a range of 
values of v for the NaRb mixture. Figure El shows a contour plot of the instability 6(q) as 
a function of Vo and the fluctuation wave-vector q. As can be seen the instability is indeed 
suppressed as v increases. In figure [3 we show a contour map of the instability 9{q) as a 
function of k and q for Vq = 1.217-Er which shows that although the k = Bloch state has 
been stabilized, the instability re-appears at finite k - a deeper optical lattice is required to 
stabilize a moving condensate. Figure |H] shows the instability 9(q) as a function of k and q 
for Vq = 2A34E R . In this case the phase separation instability is fully suppressed. The price 
for this supression is that the modulational instability at large k is much worse, with the 
dominant instability moved to higher q but with all q modes in the lowest two fluctuation 
bands being unstable once k > k/4. An alternative way of showing the same physics is to 
plot the stability boundary for k = Bloch states as a function of q and ngn{— ng 2 2) for 
fixed ngn = 0.5Er for a range of values of v as shown in figure H3 Increasing v shrinks the 
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FIG. 6: A density map of the instability 9(q) for k = as a function of fluctuation wave- vector q 
and optical lattice depth Vq (colour online). 
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FIG. 7: Instability 6{q) of a Rb-Na condensate as a function of k and g for Vq = 1.217Er (colour 
online). 

region of instability below the line as one would expect. 

Hence we have shown that two- component Bose condensates in Id optical lattices can 
exhibit both the modulational instability at large Bloch wavevectors associated with single 
component systems and the instability to phase separation at small k that can occur in the 
absence of a lattice. Turning on the lattice enhances the latter instability but suppresses 
the former. Hence such optical lattices could be used to hold immiscible two-component 
systems in intimate contact and to move them around within a trap. 
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FIG. 8: Instability 0(q) of a Rb-Na condensate as a function of k and q for Vo = 2A34Er (colour 
online) . 
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FIG. 9: The stability limit as a function for the largest unstable mode at k = for a condensate 
with fixed ngyi = 0-5Er as a function of ngu = ng<xi for values of Vo between and 7Er (colour 
online). 
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